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EQUILIBRATION METHOD FOR HIGH RESOLUTION IMAGING OF LUNG 



COMPLIANCE AND DISTRIBUTION OF FUNCTIONAL RESIDUAL CAPACITY 



5 FIELD OF THE INVENTION 

The present invention provides novel equilibration methods for imaging lung 
compliance and distribution of functional residual capacity (PRC) in the lung using 
hyperpolarized helium-3 (H He). 

10 REFERENCE TO RELATED APPLICATIONS 

This application claims priority to 60/267,282 filed February 8, 2001, herein 
incorporated in its entirety. 

GOVERNMENT INTERESTS 

15 This invention was supported in part by Grant Nos. K23 HL04486, RR02305, and ROl- 

HL-64741 from the U.S. National Institutes of Health. The Government may have certain 
rights in this invention. 

BACKGROUND OF THE INVENTION 

20 Standard methods for measuring FRC and compUance for the whole lung have been 

known for years (Meneely et al. Am. J. Med. 28:824-831 (1960); Darhng et al., J. Clin. Invest. 
19:609-618 (1940); Bums et al. Am. Rev. Respir. Dis. 130:580-583 (1984); Dubois et al, J. 
Clin. Invest. 35:322-326 (1956)). The use of hyperpolarized noble gases, such as ^He, has been 
demonstrated to be useful in the imaging of gas distribution (ventilation) in the human lung 

25 (Middleton et al, Mag. Res. Med. 33:271-275 (1995); MacFall et al, Rad 200:553-558 (1996); 
Kauczor et al, i?aJ 201:564-568 (1996);.Kauczor et al, J. Mag. Res. Imag. 7:538-543 (1997); 
Roberts et al, Mag. Res. Med 44(3):379-382 (2000); Black et al, Rad 199(3):867-870 (1996); 
de Lange et al, i?arf210(3):851-857 (1999); Altes et al, J. Mag. Res. Imag. 13(3):378-384 
(2001); Salerno et al, Mag. Res. Med. 46:667-677 (2001); Gierada et al, NMR Biomed. 

30 13(4):176-181 (2000)). However, what has not been available until the present invention, has 
been a method for high resolution imaging of the gas spaces in the lung, with vwy high contrast 
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between the signal intensity for the gas phase compared to the tissue phase, which would 
permit the combined imaging of FRC distribution and lung compliance 

The only known work in this field is related to imaging only the distribution of FRC by 
a single group Gattinoni et al. Am JRespir Crit Care Med 151:1807-1814 (1995)). The 
Gattinoni group used high resolution computer tomography (HRCT) methods to assess the 
gravitational distribution of FRC in patients with adult respiratory distress syndrome (ARDS). 
However, that work was limited by the limitations of the HRCT method. The "images" they 
obtained for distribution of FRC were fairly low in resolution, and more importantly, those 
images were only 2-dimensional. Moreover, the HRCT technique was not, and is not, capable 
of rapid measurements of FRC for the whole lung. 

Thus, a need has remained in the art until the present invention for a method providing 
rapid imaging of FRC distribution for the whole lung; true 3 dimensional imaging; and very 
high resolution imaging of lung volumes and pulmonary compliance. 

SUMMARY OF THE INVENTION 

The present invention offers novel equilibration methods for imaging lung compliance 
and distribution of fimctional residual capacity (FRC) in the lung using hyperpolarized hehum- 

3 3 • 

3 ( He) gas (H He). The imaging of puhnonary compliance combined with measuring the 
distribution of FRC is accomplished in the present invention by, for the first time, adapting and 
extending helium dilution principles and standard pulmonary function tests established for 
imaging the whole lung, to the imaging of lung volumes and compliance. The use of 
hyperpolarized ^He satisfies both of these requirements, and the present method has several 
advantages over the prior art. These include, but are not limited to, rapid imaging of FRC 
distribution for the whole lung; true 3 dimensional imaging; and very high resolution in the 
resulting image. 

The present invention provides equilibration methods for very high resolution, three- 
dimensional imaging of pulmonary compliance, and also methods for very high resolution, 
three-dimensional imaging of distribution of functional residual capacity (FRC) in the lung, 
wherein in each case, the method comprises delivering a predetermined volume of 
hyperpolarized noble gas into the conducting airways in each ventilated region of the 
pulmonary system and collecting local magnetic resonance image data therefi'om. 
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Also provided are equilibration methods for very high resolution, three-dimensional 
imaging of both pulmonary compliance and distribution of functional residual capacity (PRC) 
in the lung, wherein the method comprises delivering a predetermined voliune of 
hyperpoiarized noble gas into the conducting airways in each ventilated region of the 
5 pulmonary system and collecting local magnetic resonance image data therefrom. 

In these methods provided in the present invention, the hyperpoiarized noble gas is 
preferably hyperpoiarized helium-3 gas (H^He). Using inhaled hyperpoiarized (H^He) gas, 
three-dimensional MRI images are rapidly obtained. The H^He acts as a contrast agent, thereby 
providing unparalleled, high resolution images defining the airspaces. Because time of 
10 repetition (TR) is irrelevant to polarized gas imaging, a real-time series of images is obtainable. 

The methods further comprise dividing the lung images into as many distinct voxels as 
imaging resolution permits. Local lung volume is calculated by dividing average signal 
U intensity in each voxel by tracheal signal intensity. The average concentration of H^He in each 
ti voxel is determined by dividing an amount of tfHe in each voxel by volume of the voxel, and 
J 1 5 calculating amount of H^He in each the voxel by multiplying concentration of H^He in the gas 
' ^ space of the voxel by volume of gas space in the voxel Local PRC is determined by dividing 

W signal intensity in the voxel by tracheal signal intensity, and then multiplying by volume of the 
voxel. Finally compUance of the voxel is determined by calculating pressure difference 
between alveolar gas inside the voxel and pleural space outside of the lung, and calculating 
£3 20 local compliance in the voxel by dividing the volume of the voxel by transmural pressure 
gradient (tracheal minus esophageal pressures). Local PRC is computed by dividing signal 
intensity in the voxel by tracheal signal intensity, and then multiplying by volume of the voxel. 
The whole lung PRC is determined by summing all of the local PRC values in each voxel. 

Further provided in the present invention are the very high resolution, three-dimensional 
25 images of pulmonary compliance, distribution of functional residual capacity (FRC) in the 

lung, or a combination thereof, produced by the foregoing methods using hyperpoiarized H^He 
gas. 

In addition, there are provided systems for producing the very high resolution, three- 
dimensional images of pulmonary compliance and distribution of functional residual capacity 
30 (FRC) in the lung in accordance with the foregoing methods, comprising: means for collecting 
and processing magnetic resonance imaging data from the hyperpoiarized H^He gas-infused 
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lung, conducting airways and ventilated regions of the lung; means for dividing the lung 
images into as many distinct voxels as imaging resolution permits; means for calculating local 
lung volume; means for calculating average amount and concentration of H^He in each voxel; 
and means for calculating local FRC aad means for calculating compliance of the voxel and 
5 local compliance in the voxel. 

Additional objects, advantages and novel features of the invention will be set forth in 
part in the description, examples and figures which follow, and in part will become apparent to 
those skilled in the art on examination of the following, or may be learned by practice of the 
invention. 

10 

DESCRIPTION OF THE DRAWINGS 

The foregoing summary, as well as the following detailed description of the invention, 
will be better understood when read in conjunction with the appended drawings. For the 
puipose of illustrating the invention, there are shown in the drawings, certain embodiment(s), 
15 which are presently preferred. It should be understood, however, that the invention is not 
limited to the precise arrangements and instrumentahties shown. 

FIGs. lA and IB show a coronal sUce of polarized H^He MR images of the lung of a 
normal pig (FIG. lA) and the corresponding slice of lung volume segmentation (FIG. IB). 

FIG. 2 graphically depicts the correlation between lung volumes measured using the 
20 proposed method and those using "^He dilution technique. X axis = lung volume (in ml), 

calculated using the method of the present invention. Y axis = same measurements using "^He 
dilution technique. Based on linear coefficient, the R^ = 0.88 value shows acceptable 
agreement. 

25 DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Preferred embodiments of the present invention provides novel equilibration methods 
for imaging lung comphance and distribution of functional residual capacity (FRC) using 
hyperpolarized heUum-3 (H^He), which uses helium dilution principles to measure whole lung 
FRC based on well estabUshed standard puhnonary function tests. The principles of measuring 
30 whole lung volume with He are well established. However, H^He-based imaging is 

mtrinsically faster and gives much higher resolution in 3-D than other prior methods, such as 
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computer tomography. MRI measures the local H He concentrations at high resolution, and 
because H^He behaves physiologically just like regular He, the local measurements of H^He 
can be used to do many of the same things on a regional level that have been done with regular 
He for the whole lung, such as measuring volimie. Consequently, the use of local measurement 
5 of H^He to measure local limg volume produces entirely novel information. 

In practice, the lung is first inflated to the desired airway pressure and the endotracheal 
tube is then sealed, or in the case of human patients, the breath is held to fix the lung volmne. 
Next, a known volume of H He is introduced into the sealed or fixed lung volume and the 
H He IS mixed with the volume of air in the sealed lung. After adequate mixing, the lung is 
10 imaged with a homogeneous coil 

The lung images are then divided into as many distinct voxels as the imaging resolution 
permits, and the starting point for calculating local lung volume is the average signal intensity 
1^^;. in each voxel divided by the signal intensity in the trachea. Neglected fixmi the calculation are: 
^ (1) any signal decay due to oxygen or wall interactions; (2) any loss of helium fix)m the gas 
%^ 15 spaces due to absorption by tissue or blood; and (3) any contribution to the signal from the 
' p tissue space within the voxel. The average concentration of H He in each voxel is calculated in 
r. J terms of the amount of H^He in each voxel divided by the volume of the voxel. The amoimt of 
H He in the voxel is simply calculated in terms of the concentration of H He in the gas space of 
b] the voxel multipUed by the volume of the gas space in the voxel. Assuming complete 
u 20 equilibration of the H^He throughout the lung, the concentration of H^He in the gas space of the 
f-| voxel is exactly equal to the concentration of tfHe in the trachea. Consequently, the average 
^ ^ concentration of tfHe in each voxel is provided by calculating: 

Cone of H^He = (Cone of H^He in the trachea)x(voxel gas volume) (Formula 1) 
25 voxel volimie 

Further, assuming that signal intensity is directly proportional to the concentration of H^He, the 
gas volume (Vg) of each voxel is determined by calculating: 

30 Vg = (Svox / Strach) ^ Vvox (Formula 2) 

wherein Svox and Strach are signal intensities in the voxel and trachea, respectively, and Vvox is 
the volume of the voxel. 
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At the end of expiration (or at a tracheal pressure of zero), local lung volume is actually 
also the local FRC. Therefore, local FRC is quite simply calculated as the ratio of signal 
intensity in the voxel divided by signal intensity in the trachea, multiplied by the volume of the 
voxel (Formula 2). 

5 To calculate the compliance of the voxel, the pressure difference between the alveolar 

gas inside the voxel and pleural space outside of the lung must be determined. In breath- 
holding conditions, the alveolar pressure must equal the tracheal pressure, which is easily 
measured. Also during breath-hold, the pleural pressure will be nearly uniform (with a sUght 
correction for the gravitational field), and can be estimated by measuring the esophageal 

10 pressure. The local compliance in the voxel is then simply calculated from the volume of the 
voxel (Vvox) divided by the transmural pressure gradient (tracheal minus esophageal pressure). 

While small amounts of shunting (e.g., to the bronchial vessels) and dead-space (e.g., in 
the nasopharyx) are normal, a significant number of pathological pulmonary conditions are 
characterized by such unventilated areas. For example, blocked or impaired blood flow to an 

15 alveolus because of injury or disease, results in wasted ventilation or dead-space. However, 
because the volume of inspired H^He is larger than dead-space volume in the lung and related 
tissues, the H^He fills all of the conducting airways. In addition, a small volume of the H^He is 
delivered to each ventilated region of the lung. 

The subject of the present method is a mammal, either breathing or not. Preferably, the 

20 manraiaUan subject is human. The lung(s) of such subject may be either normal for mapping 
purposes, or injured or diseased for therapeutic or diagnostic purposes. Both halves of the 
subject's lung need not be in the same condition, and may be compared against each other. 
Altematively, a normal lung may also be used for comparison purposes. 

MRI imaging with H'^He is ideal for the three requirements of the preferred embodiment 

25 of the methods of the present invention - high resolution, high signal intensity, and rapid time 
resolution. Commercially available MRI systems, such as those manufactured by GE and 
Siemens, and improvements thereof are well known in the art, and may be used and adapted for 
these methods. 
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EXAMPLES 

The present invention is further described in the following examples. These examples 
are not to be construed as limiting the scope of the appended claims. 

5 Example L Images of PRC distribution and lung compUance in animal models . 

Confirming the feasibility and effectiveness of the disclosed method, preliminary data 
were collected using high-resolution images of H He distribution in the nitrogen inflated lungs 
of sacrificed pigs. 

After sedation, Yorkshire pigs (25-30 kg; « = 5) were intubated, md mechanical 

10 ventilation under isoflurane general anesthesia was maintained. Invasive lines were placed, and 
the animal was transported to the MRI scaimer. For imaging airspaces, H^He was produced by 
the optical pumping spin-exchange technique. An efficient and compact gradient-echo pulse 
sequence with the following imaging parameters was used: TE/TR = 0.8/2.8 msec, matrix = 
128 X 128 X 10. The H^He ventilation images of the pig lung appeared homogeneous and 

15 uiuform throughout the lungs (FIGs, 1 A and IB), making the image processing task much 
easier, and hence, more robust. 

The lung volume was measured by segmenting the lung using scale-based fuzzy 
connectedness (see, e.g., Udupa et al. Graphical Models and Image Processing, 58:246-261 
(1996); Saha etal^ Computer Vision and Image Understanding, 11:\ASA1A (2000)). Scale- 

20 based fiizzy connectedness is a region-growing technique starting firom seed points specified in 
the core of the lung. Scale is described as a local morphometric parameter, and is defined at an 
image point p as the radius of the largest hyperball inside a homogenous region with center at 
p. A hard segmentation of the fiizzy lung image was obtained by using a optimal thresh- 
holding method (Udupa et aL^ IEEE Transactions on Pattern Analysis and Machine 

25 Intelligence, 23:689-706 (2001)) that selected the threshold that best compUed with object 
morphology. 

Scale-based fiizzy segmentation of the ventilation images generated an image that 
assigned high membership values to lung tissue regions, and low values to non lung regions. 
The optimum threshold selection method generated hard segmentations firom the fiizzy lung 
30 images that were visually acceptable for all cases. The lung volume was calcxilated from each 
segmentation. 
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To validate the results produced by this method, He dilution techniques were used aad 
the total volume of the lung was measured for each case. The lung volumes obtained using the 
two methods were compared, as shown in FIG. 2. As shown, the x axis depicts the lung 
volume (in ml), calculated using the method of the present invention; the y axis depicts the 
5 same measurements using previously validated ^He dilution techniques. A linear correlation 
co-efficient was computed for this data, and the = 0.88 value shows an acceptable 
agreement. 

When summed, the measurements of regional lung volumes equal standard 
measurements of whole lung volimies. As a result, this calculation validates the techniques 

10 used in the present example against a known benchmark. This was necessary because there is 
no established technique to vaUdate measurements of regional limg volumes. (There are only 
techniques known for measuring whole lung volumes). 

To further reduce this technique to practice (1) the images of H^He distribution are 
converted into images of FRC distribution and images of lung compliance, and the sum of local 

15 FRC measurements are compared to whole lung FRC measured by standard "^He dilution; (2) 
the same experiments are repeated in the lungs of hve animals in the presence of oxygen, and 
corrected for signal decay in oxygen to show that the corrected images of FRC distribution and 
compliance matched those obtained in the lungs of a comparable pig without oxygen; (3) for 
comparison a unilateral injury will be inflicted on one lung of the test animal, thereby reducing 

20 compliance in that limg in a predictable way, then it can be determined whether the measured 
local changes in compHance match the predicted changes. Although correction of O2 decay is 
not a standard practice in imaging FRC distribution and compliance, it was necessary in 
applying the methods of the present invention because in live subjects lung compliance changes 
(decreases) in connection with various injury states, such as ARDS and pneumothorax. 

25 Because there have been no prior reported methods for imaging FRC or compHance, 

there are no existing standards available for direct comparison to the equilibration method. 
However, whole lung FRC is readily measured by standard helium dilution techniques, and the 
sum of all the individual FRC in each voxel is equal to that of the whole lung FRC, providing a 
basis for comparison with prior methods. 
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Example 2 . Single breath/pause method for measuring distribution of tidal volume and for 
imaging lung compliance . 

The single breath/pause method is also a technique for imaging limg compliance, and 
for measuring distributions of tidal volume throughout the lungs. In all subjects there were no 
apparent adverse effects from inhalation of the hyperpolarized helium gas. Blood pressure and 
oxygen saturation was continuously measured in the patients throughout the experiments. 

In the single breath/pause method, a single breath of He is delivered in early 
inspiration in a volume that exceeds dead-space volume of the lung, and then inspiration is 
stopped with a breath-hold. Because the volume of inspired H"^He is larger than dead-space 
volume, the H^He fills all of the conducting airways. In addition, a small volume of the H^He 
is deUvered to each ventilated region of the lung. Images of the H^He are obtained 
immediately while the breath is held. 

Magnetic resonance data were acquired at the MR imaging (MRI) center of the Hospital 
of the University of Pennsylvania using a superconducting 1.5 Tesla imaging system equipped 
with high-powered gradients (Echospeed Signa, 5.6 revision software, GE Medical Systems, 
Milwaukee, WI). Proton imaging was performed at 63.8 MHz using a standard birdcage body 
coil for both signal transmission and reception. Helium imaging was performed at 48.65 MHz 
using a custom-built, broadband 26 cm diameter octagonal transmit/receive surface coil placed 
on the anterior thorax. 

Pulmonary ventilation imaging was performed using hyperpolarized ^He (H^He). The 
hyperpolarized gas was prepared through spin-exchange collisions with optically-pumped laser 
rubidium atoms. (Colegrove et al, Phys. Rev, 132 :2561-2572 (1963) ; Gentile et al, Phys. 
Rev, A, 47 :456~467 (1993) ; Happer et al, Phys. Rev. A, 29 : 3092-3 110 (1984) ; Ebert et al. 
Lancet 347 :1297-1299 (1996) ; Saam, Nat Med, 2 :358-359 (1996)). Polarization levels of 
5%-15% were achieved. After polarization, the gas was transferred to the animal subject while 
the patient/subject was lying in a supine position in the MRI magnet device. The subjects 
inhaled approximately 15 cc/kg of the H^He gas to near total lung capacity (TLC). 

To ensure that sufficient magnetization would be left after all radiofrequency (RF) 
pulses used in the imaging pulse sequence, a nominal constant flip angle of 12"" was applied. 
The repetition time (TR) was selected to be as short as possible as to reduce depolarization due 
to the presence of O2. The echo-time (TE) was also chosen to be as short as possible in order to 



438312_1.DOC 



9 



CLIENT REF.N-2376 



Attorney Docket 22253-68187 Express Mail Label No. EL89 1689690US 

minimize unwanted echo attenuation due to the high diffusion of the noble gas. The imaging 
parameters were as follows: TR = 18 minutes, TE = 3 minutes, matrix size 256 x 128, field-of- 
view (FOV) = 40 X 40 cm, slice thickness = 7 mm, and number of sUces = 6. The voxel 
resolution was 1.5 x 3 x 7 mm. 
5 As above, the lung is divided into as many distinct voxels as the imaging resolution 

permits. In each voxel, the amount of tfHe that enters the voxel is diluted into the gas aheady 
present in that region of the lung. As described above, the average concentration of H He m 
the voxel is the total amount of H^He in the voxel, divided by the voxel volume Vvox- The total 
amount of H^He in the voxel is provided by the concentration of H^He in the inspired gas, Cg, 
10 multiplied by the volume of the tfHe that entered the voxel, Vg. Again, assuming that signal 
intensity is directly proportional to the concentration of H^He, the volume of the inspired H^He 
that entered the voxel is represented by Formula 2 (Vg = (Syox / Strach) ^ Vvox). 

A very simple expression relates the amount of tidal volume entering each voxel to the 
^ voxel volume multiplied by the ratio of voxel signal intensity to signal intensity in the trachea. 
^ 15 To calculate local compliance from this distribution of tidal volume, required determining the 
■ P pressure differential between the inside and the outside of the lung, which was equal to tracheal 
pressure minus esophageal pressure. This is because at the end of the breath-hold, alveolar 
pressure must equal tracheal pressure, and pleural pressure must be nearly uniform. 

Very high-resolution images of the gas spaces, with very high contrast between the gas 
C3 20 phase signal and the tissue signal are produced of the inflated lung to permit the greatest 
n possible number of voxels to be acquired. This method permits images to be acquired with 
= ^- very rapid time resolution, which is essential since the images must be acquired immediately 
after initiation of the breath-hold to avoid diffusive mixing of the H^He in the airways. 
Although presently generated in the absence of oxygen, the volume measurements have been 
25 validated. 
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Example 3 . The dynamic single breath method for imaging local ventilation . 

This method is an extension of the single breath/pause method, in which rapid, 
sequential images of image intensity are obtained in each voxel throughout inspiration as a 
5 single breath of H^He is inhaled. Each image in the sequence is analyzed according to the 
signal intensity divided by the product of the signal intensity in the trachea multipUed by the 
voxel volume. This calculation provides the volume of inhaled gas entering each voxel as a 
function of time, and provides a high-resolution image of local ventilation. 

In sum, it is feasible to obtain in vivo total lung capacity measurements and regional 
10 lung volxmies by H^He MRI thereby providing a functional lung imaging method, which adds 
potential quantitative information for assessment and management of patients with lung 
diseases. 

Each and every patent, patent application and publication that is cited in the foregoing 
specification is herein incorporated by reference in its entirety. 

15 While the foregoing specification has been described with regard to certain preferred 

embodiments, and many details have been set forth for the purpose of illustration, it will be 
apparent to those skilled in the art without departing from the spirit and scope of the invention, 
that the invention may be subject to various modifications and additional embodiments, and 
that certain of the details described herein can be varied considerably without departing from 

20 the basic principles of the invention. Such modifications and additional embodiments are also 
intended to fall within the scope of the appended claims. 
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